NASA/TM — 2003-212536/REV1 AIAA-2003-4934 



Control of Thermo- Acoustics Instabilities: 
The Multi-Scale Extended Kalman Approach 

Dzu K. Le, John C. DeLaat, and Clarence T. Chang 
Glenn Research Center, Cleveland, Ohio 


a 


October 2003 


The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key part 
in helping NASA maintain this important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the Lead Center for 
NASA's scientific and technical information. The 
NASA STI Program Office provides access to the 
NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA's institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STI Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA's counterpart of peer- 
reviewed formal professional papers but 

has less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to NASA's 
mission. 

Specialized services that complement the STI 
Program Office's diverse offerings include 
creating custom thesauri, building customized 
databases, organizing and publishing research 
results . . . even providing videos. 

For more information about the NASA STI 
Program Office, see the following: 

• Access the NASA STI Program Home Page 
at http:llwww.sti.nasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access 
Help Desk at 301-621-0134 

• Telephone the NASA Access Help Desk at 
301-621-0390 

• Write to: 

NASA Access Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076 


NASA/TM — 2003-212536/REV1 


AIAA-2003-4934 



Control of Thermo- Acoustics Instabilities: 
The Multi-Scale Extended Kalman Approach 

Dzu K. Le, John C. DeLaat, and Clarence T. Chang 
Glenn Research Center, Cleveland, Ohio 


Prepared for the 

39th Joint Propulsion Conference and Exhibit 
cosponsored by the AIAA, ASME, SAE, and ASEE 
Huntsville, Alabama, July 20-23, 2003 


National Aeronautics and 
Space Administration 


Glenn Research Center 


October 2003 


Acknowledgments 


The authors sincerely thank Dr. Daniel Paxson and Mr. Daniel Vmak for their efforts in model and simulation 
developments to validate this control design. We also greatly appreciate the contributions of the United Technology 
Research Center team to this research, especially Dr. Andrzej Banaszuk, for his invaluable insights and modeling 
data. Funding support for this work has been provided by the NASA Smart Efficient Components Project through 
the help of Project manager, Robert Corrigan, with partial support for advanced control developments by the 
Strategic Research Fund (SRF) and the Aerospace Flywheel Technology Program. 


Document Change History 


This printing, numbered as NASA/TM — 2003-2 12536/RE VI, October 2003, replaces the previous version, 
NASA/TM — 2003-212536, August 2003. The reference numbers on page 3, paragraph 4, line 6, have been 
changed. 


The Propulsion and Power Program at 
NASA Glenn Research Center sponsored this work. 


Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 


Available electronically at http: / /gltrs.grc.nasa.gov 


CONTROL OF THERMO-ACOUSTICS INSTABILITIES: 
THE MULTI-SCALE EXTENDED KALMAN APPROACH 


Dzu K. Le, John C. DeLaat, and Clarence T. Chang 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 


ABSTRACT 

“Multi-Scale Extended Kalman” (MSEK) is a novel model- 
based control approach recently found to be effective for 
suppressing combustion instabilities in gas turbines. A control 
law formulated in this approach for fuel modulation 
demonstrated steady suppression of a high-frequency 
combustion instability (>500Hz) in a liquid-fuel combustion 
test rig under engine-realistic conditions. To make-up for 
severe transport-delays on control effect, the MSEK controller 
combines a “wavelefMike “Multi-Scale” analysis and an 
“Extended Kalman Observer” to predict the thermo-acoustic 
states of combustion pressure perturbations. The commanded 
fuel modulation is composed of a damper action based on the 
predicted states, and a tones suppression action based on the 
Multi-Scale estimation of thermal excitations and other 
transient disturbances. The controller performs automatic 
adjustments of the gain and phase of these actions to minimize 
the Time-Scale Averaged Variances of the pressures inside the 
combustion zone and upstream of the injector. The successful 
demonstration of Active Combustion Control with this MSEK 
controller completed an important NASA milestone for the 
current research in advanced combustion technologies. 

1. INTRODUCTION 

Combustor instability control is an enabling technology for 
superior engine performance. Acoustic resonances driven by 
heat release fluctuations can produce large pressure 
oscillations inside the combustion chamber. 1 These thermo- 
acoustic instabilities can potentially lead to premature 
mechanical failures. In the quest to reduce particulate and 
NOx emission, future aero-engine combustors will likely have 
lean-burning front-ends. As the ground-based gas-turbine 
field has experienced, however, lean-burning combustors are 
more susceptible to combustion instabilities due to the more 
homogeneous transmission medium (an environment 
favorable for the forming of distinct acoustics modes); less 
damping from reduced or eliminated liner film cooling; and 
more energetic perturbation sources from the more vigorous 
fuel-air mixing. 

Active control of combustion instabilities by properly-phased 
fuel modulation is a very promising method for instability 
control in aero-engines. Conceptually, control of combustion 
instability pressure can be achieved by other means 2 ’ ’ 14 . 
However, these other means (e.g., speaker actuation) tend to 


be bulky, have limited operating range, or have other 
undesirable side effects. In comparison, liquid fuel 
modulation is compact, more practical, and can deliver the 
most energy with the least effort. 

The NASA Glenn Research Center (GRC) Active Combustion 
Control (ACC) technology effort aims to demonstrate active 
control in a realistic and relevant environment by providing 
experiments tied to aircraft gas-turbine engines. The first 
successful demonstration of active control of a High 
Frequency (HF; >500Hz) thermo-acoustic instability via fuel 
modulation under engine-realistic conditions was 
accomplished in June 2002 by the NASA GRC effort. It was 
reaffirmed in September 2002 with additional experiments and 
improved results. These experiments were done in 
collaboration with Pratt Whitney 14,15,16 at the United 
Technology Research Center (UTRC), on the NASA Single 
Nozzle Combustion Rig (SNR). 

Prior to this work, partial successes with suppression via fuel 
modulation of Low Frequency instability (LF; -280 Hz) had 
been reported elsewhere, including the works supported by 
NASA. 26,27 In these earlier attempts to suppress the dominant 
instabilities, other instability peaks would arise around the 
central frequency (so-called “peak splitting” 26,27 as seen on the 
pressure amplitude spectra plot, Figure 1). A solution to this 
problem, which occurs even with moderate control gains, has 

not been identified in previous works 8 14 . For this effort, 

active instability suppression has been demonstrated for 
high-frequency oscillations on the NASA SNR. This HF 


Combustor LF Instability - Run 1 14, Pt 55 and 42, 8/19/2000 



Frequency, (Hz) 

Figure 1: Low-frequency instability suppression 
control, showing peak-splitting 26 
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Figure 2: Engine and NASA Single Nozzle Combustion 
Rig Comparison 

instability in the NASA SNR represents an actual case of 
engine instability although the rig exhibits higher background 
noise/disturbance levels than the engine (Figure 2). The 
NASA SNR operates at engine pressures, temperatures, and 
flows using liquid jet fuel and incorporates many engine-like 
features including an actual engine fuel injector (Figures 3 and 

4 ^ 14 , 15,16 

This successful demonstration of ACC for the HF case was 
carried out with two distinct controllers developed by NASA: 
an adaptive controller by Kopasakis 17 and an observer-based 
controller by Le in the Multi-Scale Extended Kalman (MSEK) 
approach described herein. The two controllers yielded 
comparable results, and neither showed any sign of “peak- 
splitting”. This paper is devoted to the MSEK model-based 
controller. 

In this paper, a brief introduction to the MSEK control 
approach is given. Next, the MSEK control design approach 
is illustrated with a detailed description of the controller 
developed for the experimental demonstration of combustion 



Figure 3: Combustor Rig for Combustor Instability 
Control Research (Photo by United Technologies Research Center) 
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instability control. Results of the evaluation of the controller 
via numerical simulation are presented followed by the 
experimental results with the combustor rig. Finally, a 
generalization of this MSEK controller and some concluding 
remarks are given. 

2. THE MULTI-SCALE EXTENDED KALMAN 
CONTROL APPROACH 

State variables tracking and optimal control normally can be 
handled effectively by Kalman Filters and the LQR method 18 
only if the system dynamics are fairly close to the classical 
case of stochastic processes. These classical control methods 
are often found inadequate for aero-engines or other complex 
systems with a high level of dynamic uncertainties, complex 
(external or internal) disturbances, noise, and transport delays 
in the control mechanisms. Examples of such complex 
disturbances include: mass-imbalance; sensor distortion; high- 
order or nonlinear behaviors of rotor, motor, and bearing 
dynamics; unknown excitations of thermo-acoustics 
instabilities; and higher-order stall and surge dynamics; etc. A 
feasible solution for many such difficult cases is to predict the 
variations in system transient characteristics with Multi-scale 
Analysis (e.g., with the Wavelet Transform 19 ’ 20 ’ 21 ’ 22 ) to adjust 
model and control parameters accordingly, and to compensate 
the control commands for unknown, transient disturbances. 

That direction for control design was proposed by Le 21 and 
named the Multi-Scale Extended Kalman approach in this 
paper. It was inspired by the fact that, non-stationary 
processes (e.g., the Fractional Brownian Motions) become 
stationary under proper wavelet transformations. 20 So, Multi- 
Scale Analysis promises a relief from the limitations of 
Kalman filters, which provide optimal estimation solutions 
only to stationary processes (e.g.; linear system driven by 
white-noise — a stationary random process). The mathematics 
foundation for Wavelets representation of (non- trivial) random 
processes was established by Houdre 20 in 1993. Implicit in 
Houdre results is the advantage of minimizing estimation 
errors or system outputs variance after a proper wavelet 
transformation of the related random processes onto the time- 
scale domain. Thus, the ensemble variance (or a time- 
averaged approximation of the variance) of system outputs -- 
the usual optimization criteria in the classical Kalman 
approach — would be replaced in the MSEK approach by the 
time-scale averaged variance (limited to selected ranges of 
signal details or “scales”) of the transformed signals. 

The MSEK method has been found to be effective for many 
difficult problems where it is necessary to track system states 
with high accuracy despite disturbances, system uncertainties, 
and non-linearity. This approach requires state-variable 
estimation for the system dynamics specially formulated in the 
general framework of an Extended Kalman Observer - i.e., 
extensions of the Classical Kalman Filter technique to handle 
the nonlinear and time-varying natures of complex systems 23 . 
In addition, wavelet-like, Multi-Scales models instead of the 
ordinary Fourier bases are integrated into the Extended 
Kalman Observer to estimate system disturbances. The 
inclusion of these Multi-Scales disturbance models and Multi- 
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Figure 4: Schematic of the NASA SNR Combustion Rig installed at UTRC 


Scale analysis into the Extended Kalman system state 
equations allows for effective suppression of transient 
disturbances as well as adjustment of observation and other 
control parameters against system uncertainties. 

Multi-Scale signal-processing techniques are versatile and 
powerful for analyzing both transient and steady dynamics. 
Wavelet bases, for example, provide coherent building blocks 
for signal analysis and synthesis. They, however, are too 
computationally intensive for use in real-time monitoring or 
control. Approximations and simplifications of various 
wavelets and associated algorithms are necessary to reduce 
computation time. Nevertheless, for many challenging 
applications, such customized implementation of Multi-Scale 
techniques still provides greater versatility and accuracy than 
does the Fourier Transform approach. Practical wavelets such 
as the Compact Harmonic Wavelets 22 are effective for 
estimating and compensating for quasi-periodic, complex 
disturbances, or for fine-tuning Kalman filters. 

So, in this method the Multi-Scale digital filters are used to 
calculate the “Time-Scale- Averaged Variances” of system 
outputs as performance metrics for real-time tuning of model 
and control parameters. Such Multi-Scale Performance 
Metrics effectively differentiate competing parameters (with a 
generally well-conditioned Jacobian on the parameter space), 
allowing for sequential search and adjustment of their optimal 
values. 

This new control design approach and several practical Multi- 
Scale analysis algorithms were developed for system 
identification and real-time tuning of model and control 
parameters, 21, 22 and have been successfully applied since 
1999 to the control of magnetically levitated, high-energy 
flywheels. 24,28 Several test rigs of representative problems in 
engine and rotating machinery were used in the development 
of various techniques to combine Multi-Scale Analysis with 
Extended Kalman filters to estimate transient disturbances and 
enhance system states tracking. 


3. A MSEK CON FROL DESIGN FOR ACC 

A Phase-Shifting Proportional scheme alone, even one with 
phase-drift adjustments, is only effective against relatively 
steady harmonic excitations. For combustion instability 
control, the highly transient thermal excitations and other 
unknown disturbances will induce side-lobes (“peak-splitting') 
and hence severely reduce the effectiveness of the controller in 
damping out the acoustic modes. 27 This motivates the 
following control design strategy: 

(a) A Predictive Extended Kalman Filter to reduce the 
adverse effect of large transport-delays and random 
disturbances for effective damping of the acoustic 
modes. 

(b) A parallel Tone Suppression scheme with varying 
gain and phase to counter thermal excitations and 
other unknown harmonic disturbances. 

This paper will describe the numerical validation and 
empirical demonstration of the control strategy as summarized 
by (a) and (b). A more thorough presentation of the 
mathematical fundamentals underlying the MSEK control 
design is expected to be published in the near future. 

The algorithms and implementation of (a) and (b) are 
presented in the next three sub-sections, beginning with a 
functional overview. Figure 5 shows a top level functional 
diagram of the digital controller for ACC named Multi-Scale 
Predictive Dampe r (MSPD) designed in the MSEK approach. 

3.1. Design Functional Overview 

The greatest challenge for combustion instability control is to 
overcome the inherent, large time-delays on fuel modulation 
effects under severe disturbances (random as well as 
harmonic). The total time-delay could be longer than three 
oscillation cycles of instability. Thus, it is necessary to have 
a Kalman-type predictor for the dominant acoustics states -- 
Sub-block “ EK States Predictor ” in Figure 5. In the Predictor, 
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Figure 5: Multi-Scale Predictive Damper (MSPD) functional block diagram 


the modeled dynamics of combustion pressure oscillations 
include the effects of (non-linear) heat release fluctuations. 
Inputs to the Predictor are not the combustion zone sensed 
pressure, but rather the ‘"tones” extracted from it by a Multi- 
Scale Tones Analysis scheme, as will be explained in the 
following sub-sections. The fuel modulation command 
includes a proportional component intended for damping out 
the dominant acoustic modes based on the Predictor and a 
non-linear component ( Tone-Suppression ) for countering 
thermal excitations and other transient disturbances. These 
excitations are estimated by a phase-adjusting Multi-Scale 
tones extraction scheme ( Phase- Adjusted Reconstruction ). 
The adjustments for phase drifts are calculated by a regression 
scheme ( Phase-Drift Estimation ). Besides the phase-adjusted 
tones, the Multi-Scale Analysis function in MSPD also 
calculates the Time-Scale Averaged Pressure Variance which 
is used as a performance metric for automatic global search of 
optimal parameter ranges ( Parameters Tuning ) and 
adjustments (e.g., control phases and gains, and predictor 
parameters). The compensation for undesirable pressure 
disturbances from the fuel line just upstream of the injector 
( Upstream Compensation ) is also a damper. 

3.2. Tone ‘"Quadratures” Extraction 

The purpose of the Multi-Scale Tone Analysis block in MSPD 
is to extract predictable, transient characteristics in the 
pressure measurements for control input. It was found that, 
the Single Nozzle Combustion Rig (SNR) might exhibit either 
the LF (~285Hz) or HF (>500Hz) oscillations depending on 
the length and location of one of the two spool sections 
(upstream airflow) relative to the pre-diffuser. Figure 4 shows 
the HF configuration of the Combustion Rig, for which both 
spools are upstream of the pre-diffuser that flows into the 
combustor dome. The Tone Analysis part of controller MSPD 
(Figure 5) is performed by wavelet-like filters to extract the 
wavelet ‘"quadratures” of the “tones” in the sensed pressures, 
in a vicinity of the two frequencies. 


By definition, a “quadratures” pair of a signal stream s(t) in a 
frequency band is any decomposition of the signal contents of 
that band into two orthogonal components Si(t) and s 2 (t) by 
using an Orthogonal Wavelets Decomposition and 
Reconstruction scheme (e.g., using “quadrature mirror 
filters”). 19 In this strict definition of signal “quadratures”, the 
mutual- orthogonality condition imposed on S|(t) and s 2 (t) is 
understood in the sense of L 2 -Function Inner-Product and is 
automatically satisfied even for the finite support intervals of 
the Orthogonal Wavelets of choice. 

For simplicity, instead of a formal Orthogonal Wavelet 
Decomposition and Reconstruction scheme, a pair of Finite 
Impulse Response (FIR) band-pass filters is designed for each 
frequency band to extract the signal (approximate) 
""quadratures”; one filter for extracting sj(t) with a phase lead 
of 45 degrees at the center frequency of interest, and the other 
for s 2 (t) with a phase lag of 45 degrees; the gains of both 
filters are equal to 1 at the central frequency. The advantage 
of these multi-scale (approximate) “quadrature filters” is in 
their simplicity; also, the polar angle of the vector [si(t) s 2 (t)]' 
([...]' denotes a transposed matrix) represents the 
instantaneous “phase” angle of the extracted tone plus 7i/4 
radians. 

Notes : (1) In other Multi-Scale Analysis applications, 
“Quadrature Filter Banks”, such as the Compact 
Harmonics Wavelets, 22 are needed for extracting multiple 
tones from the raw signals. 

(2) The “Wavelef’-Fourier-Analysis connection: 
Roughly speaking, a “frequency band” of the Fourier 
Transform corresponds to a “wavelet subspace” 
composed of a certain range of “feature scales” (or signal 
“details”) of the Wavelet Transform. The higher the 
frequency band, the smaller the corresponding feature 
scales are. 
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3.3. Quadratures Prediction 

Throughout this paper, S„ denotes the 4x1 vector of the 

of tone 1 and tone 2 extracted 

from the sensed combustor chamber pressure oscillations (LF 
and HF tones, respectively; with the superscript as the tone 
index and “n” as the time index). The following approximate 
“evolution model” (Equations (1-a) to (1-f)) for the 
combustion pressure quadratures at N sampling steps (~ 
control-delay time) ahead of the “current” time index “w” was 
used for the formulation of the predictor (See Appendix: 
Predicted Quadratures Evolution Model): 


‘quadrature pairs” S n = 


s'„ 

si 


S ;V+ „ +) = AS„ + „ + Q+w 

z„.,=uz„ + fmS v ,„ 
y« =Ez » +DM W v 


~P\ 


’ h ll 

0 

hi 2 

0 “ 

S + 

~b 

Pi 


h 2l 

0 

^2 2 

0 

n 

b 


with limits: - p < p: < p 

y max ri r max 


Q H [o 



(for / = 1 , 2) 


0 , 0 Q 2 ]' 


d-a) 

(1-b) 

(1-c) 


d-d) 


(1-e) 


(1-0 


The “output measurement” in these equations is a 2x1 vector 
y„ representing the combined quadratures of mode 1 and mode 
2 extracted from combustor pressure measurement at time 
index “w” for input to the Kalman predictor. That is, 


y„ = ms. 


M = 


1 0 1 OT 

0 10 1 


(2-a) 

(2-b) 


The large time-index shift from Sy+ n to the feedback y n was 
approximated by a 5th-order Pade approximation. That is, in 
Equations (1-b) and (l-c),y w was viewed as the result of a 5 th - 
order Pade-transform (for a time-delay ~ N*T; T denotes the 
sampling time) applied on the 2x1 vector M S N + n . So, U, F, 
E, D are just the combined state matrices of the pair of Pade- 
transformations used (with the respective dimensions 10 x 10, 
10x2, 2x10, and 2x2). The quantities v and w are band- 
limited random vectors. 


A = 




cos(<qT) -sin^T) 
sin(a>,T) cos(cq T) 


[ 0 ] 


[ 0 ] 


4 


COS(rt*T) 

sin(<y,T) 


-sin(tf> 2 T) 

cos^T) 


( 3 ) 


with T being the sampling time; CO\ . Ob (in radians/sec) are 
the two tone frequencies; ^ / , ^2 their open-loop damping 


coefficients; and A.j=exp (-2^ \ 0)\J ); A. 2 =exp (-2£ 2 M 2 T ). 
The effect of combustion heat fluctuations on these acoustic 
quadratures (at N steps ahead of n) is approximated by the 
vector-valued function Q of S„ and the fuel modulation 
command U n , as formulated in (1-d) to (1-f)- The parameters 
J h i ■ | are constant, internal feedback gains; b the 

approximate control authority; P a a reference pressure 
comparable to the mean pressure drop across the injector. The 
parameters { h , . j, t>, /f , /? ma x » anc ^ N are chosen to match 

the open-loop forced response characteristics of the 
combustion chamber pressure. 

The Kalman Quadratures Predictor : Equations ( 1 -a) to ( 1 -f) 

result in an estimation for S_\r+ n , denoted by S fl , formulated 
explicitly in Equations (4-a) to (4-e) below: 


^ = AS„ + O + K,( EZ„ +DMS„ - MS„ ) 

Z„,=UZ„ +FMS„ +K,( EZ„ +DMS„ - MS,, ) 


P 1 
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0 h 
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P 2 _ 
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22 0. 

n 

with limits: -p 
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n / \ 




Q: = K 

Jr 

Pi / 
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1 

It 

/Po) 
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Q = [ 

0 

Q, 

0 Q 

2P 



(4-a) 

(4-b) 


(4-c) 


(4-d) 


(4-e) 


The matrices Kj and K 2 (of dimension 4x2 and 10x2, 
respectively) are the two vertical sub-blocks of the Kalman 
matrix computed from the evolution equations (1-a) to (1-c) 
for certain covariance matrices assumed for the random 
vectors v and w in these equations, treating Q as a known 
input at time Note that, in this simplified estimator for 
Sn+„, the Kalman innovation term 18 was omitted. The 
numerical values of the parameters in Equations (4-a) to (4-e) 
will be presented in Section 4. 

3. 4. Phase Bias and Drift Correction with Quadratures 


If S =[sj(/) 5j(/)]' is the vector of quadratures extracted from 

a signal S(t) using “Quadrature Filters”, then a combination 

S a = (^(/)c°S^-^(/) s inQr) represents the tone of s(t) in the 

Quadrature bandwidth phase-shifted by “a+7c/4” radians. 
Thus, it might be used in a Proportional Phase-Shifting control 
scheme for a suitable phase bias “a”. But, a Proportional 
Phase-Shifting control with a fixed phase bias is not effective 
in handling large time-delays and phase drifts of transient or 
time-varying dynamics, especially when the Signal-to-Noise 
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ratio is small. To remedy this, a regression scheme was used 
to estimate the drifts in the accumulated phase of each tone 
over the expected time interval of the total transport-delay. 

3.5. Damper Command 

A part of the fuel-modulation command u n , essentially an 
LQR damper (distinguished by the superscript ‘cT ) of fixed 
gain, is formulated as: 

Un'=Ml ^(O+ZA ^ 2 (0’ ( 5 ) 

where 5^ (/„) and (t n ) are the signals constructed 
from each of the predicted quadratures as 2><1 sub-blocks in 

a r ci 

c _ 

0/7 = — , 

C 2 

the tone indices i = 1 , 2 ), where /? is a fixed bias, 8 1 n the 
predicted phase-drift. The quantities /u\ and are fixed gains. 
The drifts 8 . n (over the time interval of expected control 
transport-delay, for tones i = 1 , 2 ) are predicted by a 
regression scheme applied on the quadrature phases {a . ; n } as 
defined in Equation ( 6 ) of the next sub-section. 

3.6. Suppressions of Residual Harmonic Disturbances 

This part of the control scheme is to cancel out remaining 
harmonic disturbances that can not be suppressed by just the 
“damper”. For this purpose, the instantaneous amplitudes and 
phases of the combustion pressure near the LF and HF are first 

estimated by projecting the quadratures S' n (for / = 1 , 2 , as 
2 x 1 vectors) onto an orthogonal frame rotating at the 
respective speed 0> x . To reduce the effect of random noises, 
these projections are next passed through low-pass filters with 
cut-off frequencies at 1 00 Hz and 200 Hz, respectively, before 
computing the amplitudes and phases. Note that these cut-off 
frequencies are about 1/3 of the respective central frequencies. 
Consequently, the predicted pressure oscillations for the 
respective frequency bands are given by: 

cr ; = r in cos (&>. t n + a , n ) (for i- 1 , 2 ) ( 6 ) 

wherer )n is the amplitude, and {a f n } the polar angle at time 
t r 7 of the filtered, projected quadrature vector. 

Note that, (J j should be highly coherent with the unknown 

“residual” harmonic excitations other than white-noise 
disturbances and the heat quantities associated with expression 
(1-f). Therefore, to counter this unknown effect the fuel- 
modulation command u n is augmented with a compensation 
term (distinguished by the superscript ‘c’) of the following 
form, as computed in Tone Suppression : 


using respective phase angles /?. = P+8 in (for 


W n ' = X r> n r i nCOS (d), t„ + « in + (X+S J 
/' = 1 , 2 

(for tones / = 1, 2), where y jn is an varying gain applied on 

(Jj after it is phase-shifted by a bias a plusr? /n . These 
gains are calculated as discrete values of this expression: 

Yi =K + 7 jift ~P;)dt (for i=l,2) , (8) 

where y () is a constant gain compatible to the expected control 
authority; 77 , a constant integral gain; pi , a reduction level 
expected to be maintainable for r, with this scheme. Also, y^’s 
are limited between certain bounds y min and y max before 
substituted into Equation (7). 

3.7. Fuel Modulation Command and Tuning 

Finally, expressions (5) and (7) combined gives the fuel- 
modulation command: 


K = *£' + < 


(9) 


In this control scheme, all the parameters in the predictor (as 
well as the control scheme) appear to be fairly deterministic 
and empirically predictable (at least for the SNR case), except 

for the phase biases, a and P , and the control authority 
represented by “b” (see Equations (1-d)). These parameters 
are set by the Parameters Tuning logic in sequential global 
searches at low loop-gain. Once found, the average phase 
biases are adjusted with a slow- varying adjustment term 8. 
Such setting and adjustment of control variables are based on 
minimizing the Time-Scale Averaged Pressure Variance , 
defined for this controller to be a 0.2-Hz low-pass version of 
the square-norm of S n . 

Note : The pressure perturbations just upstream of the 
injector are highly correlated with the combustion 
pressure. So, they can be used as additional feedback to 
provide more active damping to the acoustic modes either 
inside the combustion zone or in the fuel-line. This can 
be done with a predictive damper as in the calculation 

of u {I . Such dynamic equations for the predictor are far 
from being a physical model of the fuel-line dynamics. 
But, the actual correlation of the fuel-line with the 
combustion zone dynamics should nevertheless be 
adequately retained in the predicted states by virtue of the 
Kalman gain matrix. 


4. DESIGN EVALUATION AND EXPERIMENT 

The NASA Single Nozzle Combustor Rig (SNR) at UTRC 
was used as the experimental test bed for demonstrating 
closed-loop instability suppression using the MSPD controller. 
The SNR was operated at a pressure, temperature, and fuel-air 
ratio corresponding to a mid-power engine condition 
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(T3=770°F, P3=200psia, fuel-air ratio=0.03) and exhibited 
roughly the instability behavior shown in Figure 2. 

4.1. Simulation Models 

A reduced-order, parametric model and a “Sectored- 1 D” 
thermo dynamics flow model (Paxson 2:> ) of the SNR were 
used in the design and evaluation of the MSPD controller. 
The former is essentially a Helmholtz-type model of the two 
thermo-acoustics modes of the SNR, formulated in quadrature 
state variables similar to Equation (1-a) to (1-0 as follow: 


T h on control, phase drifts due to some randomness in 
the delays, and disturbances W n . 

4.2. Control Validation Results 

The following parameter values of the reduced-order model 
for the HF-Configuration SNR are with PSI as the unit for the 
pressure feedback, Volt for the fuel-modulation command u n 
to the servo valve used — which was limited to 2 Volts in 
amplitude. 


X n + \ =AX n +0(T a X n ,T b U n , T c p„,W „ ); (10-a) 


p° c = [c, 0 c 2 0]X n 

S n =P°n C ' +M, 


hi i 0 hj 2 0 
h 2l 0 h 22 0 


dx,Y 


b g 1 
b g 1 


rlJ ai ) 

c n 

w 

n 


with limits : —p < n < p 

max max 


Q: —K 


Pi , 


f Po 


+1 -1 


(for i = 1,2) 


Q s [o Q, 0 Q 2 ]' 


(10-b) 

(10-c) 

(10-d) 


( 1 0-e) 


(10-0 


In equations (10-a) to (10-f), the state matrix A is as in 
expression (3). And, 


y _ r 1 l 2 2-.* 

A. n — yx x X X J stands for the 4x1 state vector 

1 n 2n \n 2n 


First, the measurement gain and the internal feedback gain for 
the dominant mode (Equation (10-b) and (10-d), respectively) 
are both normalized so that: 

C 2 2 = 1 ; and h-i 2 = \ 

With that normalization and units choice, the rest of model 
parameters are chosen to match empirical data for the 
sampling time T = 0.0002 seconds (i.e.; 5000 Hz I/O rate are 
used for both the reduced-order model and controller outputs 
on the rig): 

C, , = 0.8 ; 

#,=20 ; K=0.6p o =l2 ; 

h, | 1 ; h, 2 — h 2 1 1 > 

T a =0.0038; T b =0.0033; T c =0.0002 (mean values); 
b =0.2322 ; g =0.01 ; 

The rms value of W n is 0.2805 for the high disturbance 
setting, and is 0.07 for low disturbance; that of V n is 0.07. 


QC 

of mode 1 and 2 (as superscripts); p n the static pressure 
oscillation about the mean static pressure p (of the operating 

' ac' 

condition in the combustion zone); S fV sensed p n ; and 

W n , V n are scalar, white disturbance and sensor noise 
limited to below 2000Hz. The heat-release fluctuation effect 
Q is a vector-valued function of delayed states, fuel 

modulation command, U n , including valve dynamics, and 

disturbances. The operators T a , T b , r. are random time 

shifts (with non-zero mean values) applied on the respective 
variables. The effect of the combustion zone pressure 
oscillation on fuel modulation is too complex to model 

accurately; so, the terms involving p“" and the constant gain 

“g ” are used here only for testing the robustness of the 
controller, not to be included in the predictor equation. 

Remark : This plant model presents all the main 
challenges of the SNR: the thermal excitations modeled 

by the internal feedback matrix { h ; }, the time delay 


The effect of fuel valve on the command a n is represented by 
a low-pass filter (resulting in U n ) and an attenuation which is 
already absorbed into b. 

The damping parameters in the predictor (Equation (4-a) to 
(4-e)) and in Equations (1-a) to (1-c) for the Kalman matrix 
calculation were deliberately set at 30% of critical damping 
for both modes, not at the estimated values which is only 
about 0.03 for the HF mode. This is to ensure stability for the 
predictor with a simplified implementation of the Kalman 
feedback. 

The validation results for the MSPD controller using this 
reduced-order model of the SNR are summarized in Figures 6 
to 8. Figure 6a presents the Amplitude Spectra and 
Probability Density of combustor pressure oscillations with 
controller-on in comparison to the controller-off case, both 
with a much lower level of disturbances than in the SNR; 
however, this disturbance level is compatible to a known 
actual engine environment. Figures 6b shows the favorable 
time- domain response of combustor pressure to the controller 
for this low-disturbance case, as the controller is turned ON 
and OFF alternately every 5 seconds. Figure 7 further 
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(a): Closed-loop versus open-loop combustion pressure 
oscillation Amplitude-Spectrum and Probability Density 


(b): Pressure oscillations (top) and commands ON/OFF 
(bottom) 


Figure 6: Simulation results with the reduced-order model at low random disturbance 


demonstrates the absolute stability achieved by the controller 
(despite the large delay-time in control action) when there is 
no external disturbance, with the controller turned ON and 
OFF subsequently every 5 seconds. Lastly, Figure 8 presents 
the Amplitude Spectra and Probability Density of combustor 
pressure oscillations with controller-on in comparison to the 
controller-off case, both with a high level of disturbances 
comparable to the SNR environment. 

4.3. Experimental Results 

For the experimental demonstration of the MSPD controller, 
combustor pressure in the SNR was sensed about 2 inches 
downstream of the fuel injector. The control algorithms were 
implemented on a real-time processor. The fuel flow was 
dynamically controlled via a high-response fuel valve. 


MSEK ACC vmUMbon without random disturbance) (Control ON/OFF) 


ifj 

.l iy Muk iya 

* Wr 

■’o 0 % 

? i 

1 tt 

2 2 5 3 3 5 4 4 5 5 
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2 l 1 

| ON Off 

c A , 

ON OFI 




E v ~ 

I- 

-si 1 






0 05 1 15 2 25 3 35 « 45 5 

Tune (second*) 


Figure 7: Pressure oscillations (top) and commands 
ON/OFF (bottom); simulated with the reduced-order 
model without disturbance 


The significance of the “Predictor-based Damper” in the 
MSPD controller is mainly its effectiveness in supplying 
active damping to the acoustics modes. It is expected to be 
highly effective when the “background” disturbance level in 
the combustor is comparable to an engine environment such as 
Figure 2, not four times higher in Amplitude Spectra as in the 
SNR. But, it becomes less effective in the noisy environment 
of the SNR. This is when the "Tone Suppression" function in 
the controller becomes necessary. With Tone Suppression, the 
controller performs well, as predicted on the reduced order 
SNR model for the September 02 tests. On the other hand, 
although the nonlinear “Tone Suppression” part of MSPD has 
the most powerful effect on the PSD center peak, it needs the 
Damper to prevent “peak-splitting”. 


j03 

g02 


Amplitude Spectrum - MSEK ACC si mutation at high random dsturbances 


Open-toop 

(with high random dsturbances) 
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200 300 400 500 600 700 800 
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Figure 8: Closed-loop versus open-loop combustion 
pressure oscillation Amplitude-Spectrum and Probability 
Density; simulated with the reduced-order model at a high 
random disturbance 
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Combustor Pressure Amplitude Spectra, psi 


Figure 9 shows the Amplitude Spectra of combustion pressure 
oscillations in an experiment with the MSPD controller on the 
Single Nozzle Combustion Rig; for comparison, this result is 
overlaid on that of the baseline instability (controller turned 
off). This result agrees with the earlier prediction on the 
reduced-order model. Both show about 30% reduction of 
pressure oscillation amplitude at the HF peak when the MSPD 
is used. The Upstream Compensation feature in the controller, 
however, has not been tested on the SNR. 

Remark : Interestingly, the technique used in “Tone 
Suppression” was found also effective for mass- 
imbalance compensation in controlling magnetically 
levitated flywheels. 

5. GENERALIZED MSEK CONTROLLERS 

As mentioned previously, Houdre established the mathematics 
foundation for Wavelets representation of (non-trivial) random 
processes in 1993. 20 A practical framework suitable for 
engineering applications of Multi-Scale-transformed random 
processes was introduced by Le 21 in 1995. Continuing in this 
direction, Equations (1-a) to (1-f) lend a useful example of 
Multi-scale building-blocks for the approximate formulation 
of system dynamics in Multi-Scales variables. The advantage 
of such approximate model is in the formulation of Kalman 
predictors or estimators of system transient characteristics 
with enhanced sharpness via Multi-Scale analysis. 
Quadratures evolution equations like (1-a) to (1-f) are suitable 
for the modal dynamics of acoustics, structures, rotors, etc. In 
general, the dynamics of multiple tones/modes can be 
approximated as such, using Compact Harmonics Wavelets 
Filter Banks, 22 but probably with slightly more complex 
internal feedbacks (e.g., Equation (A4) in Appendix) and 
requiring more sophisticated. Extended Kalman filters. For 



Figure 9: Amplitude Spectra showing the effects of active 
combustion control by the MSEK method over the combustion 
instability peak pressure oscillations 
(9-11-02 SNR test) 


more deterministic systems, such as in flywheel magnetic- 
bearing control, physical-based (extended or conventional) 
Kalman filters and Multi-Scale techniques can still be 
integrated as in Figure 5 to effectively handle non-steady 
disturbances 24 . Furthermore, many types of modal dynamics, 
including those driven by non-stationary disturbances, could 
be approximated in such a way that optimal state-variables 
observers (to minimize the estimation errors covariance on the 
time-scale domain) exist in closed-form . 

6. CONCLUSIONS 

The Multi-Scale Extended Kalman (MSEK) approach has 
provided numerous capable, novel controllers for a wide range 
of control challenges in current research at NASA in 
innovative aerospace engines and power systems. The key to 
the versatility and effectiveness of the MSEK control 
approach is its innovative combination of Extended-Kalman 
Observers and Multi-Scale “features extraction” techniques 
for disturbance estimation, compensation, and performance 
optimization in real-time. In a recent demonstration of Active 
Combustion Control (ACC) with fuel modulation on a liquid- 
fuel single-nozzle combustor rig, a MSEK controller 
suppressed high-frequency (>500Hz) combustion pressure 
oscillations by 30%. The combustor rig emulates an actual 
case of engine instability and operates at engine-realistic 
conditions. This demonstration and further analytical studies 
show the feasibility of the application of the MSEK approach 
to Active Combustion Control, which is an enabling 
technology for future aero-engine combustors with improved 
efficiency and reduced emissions. 

APPENDIX 

Predicted Quadratures Evolution Model 

Each of the two possible acoustics modes in the NASA Single 
Nozzle Rig combustion chamber (~285Hz and >500Hz) is 
approximately a second order dynamics with low damping 
(about 3% of critical damping for the dominant mode) and 
driven by the fluctuating heat release. The heat release has a 
very large time-delay for the desired control bandwidth, about 
0.0035 seconds or more, as empirically determined. With 
added convection delays, the time lapse between sensed 
chamber pressure oscillations and the associated heat 
fluctuation could be more than 0.0038 seconds. Time-lag and 
saturation in heat-release are the important contributors to the 
sustained combustion acoustics. The effect of heat-release 
fluctuations on the combustion acoustics stays within certain 
determinable limits. Based on those assumptions, the 
following approximate “evolution model” for the combustion 
pressure quadratures at N sampling steps (~ control-delay 
time) ahead of the “current” time index “/7 ” is assumed for 
the formulation of a predictor: 

S v + „ + , = AS v+ „ + Q(S„ +jt ,M„) + W (A 1 -a) 

y„ = ^ (Al-b) 
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(A5) 


In Equation (A 1 -a) and (A 1 -b), S„ denotes the 4x 1 vector of 


the “quadrature pairs” S n = 


s: 


of tone 1 and tone 2 


(—285 Hz and >500Hz, respectively); the superscript is the tone 
index and “ n ” is the time index. 


( Note : For simplicity, instead of a formal Orthogonal Wavelet 
Decomposition and Reconstruction scheme, a pair of Finite 
Impulse Response (FIR) band-pass filters can be suitably 
designed for each frequency band to extract the signal 
approximate “quadratures”: One filter with a phase lead of 45 
degrees at the center frequency of interest, and the other with a 
phase lag of 45 degrees; the gains of both filters are equal to 1 
at the central frequency. The advantage of these “Quadrature 
Filters” is in their simplicity. The polar angle of the vector of 
extracted quadratures [S](t) s 2 (t)]' by such filters represents the 
instantaneous “phase” angle of a tone plus 7i/4 radians — [...]' 
denotes a transposed matrix.) 

The “output measurement” in Equation (Al-b) is a 2x 1 vector 

y representing the combined quadratures of mode 1 and 
J n 

mode 2 extracted from the sensed combustor pressure at time 
index So: 


M — 


10 10 " 
0 10 1 


(A2) 


The quantities W is a band-limited random vector. The 4x4 
“state-matrix” A in Equation (Al-a) represents only the 
dynamics of combustion chamber acoustics in modal 
representation, not including thermal excitations. Explicitly, 


4 


A = 


0 


0 


* 


cos(a| T) -sinfc^T) 

sin(<^T) cos(qT) L ' (A3) 

cos(<y,T) -sin(&>,T) 
sinfeT) cosfeT) 
with T = 0.0002 seconds, the sampling time; and, 

0 )\ = 570 n radians/sec; C , , = 0.30; exp (-2 £ , <y,T ). 

Oh> 1000 k radians/sec; £2 = 0.03; 1 = exp (-2^ 2 0hJ ). 


The effect of combustion heat fluctuations on these acoustic 
quadratures (at N steps ahead of n) is modeled as a certain 
vector-valued function Q of the earlier acoustic states denoted 
by S k+ n and the fuel modulation command u n (with k being a 
fixed integer much smaller than N). This vector-valued 
function is assumed of the form: 

Q( S k+„>«n) = [0 Qi 0 0 2 ]', (A4) 

where the time-varying entries on the rhs are calculated using 
the following expressions of “effective modal excitations” 

/?, and p 2 on fuel-flow variations and the acoustic states at 
time ‘W+«” : 


P\ 


"h lt 

0 

hi 2 

0" 

Pi 


h 2 i 

0 

h 22 

0 


S n+k + 


where { h ■ | are constant, internal feedback gains, b the 

approximate control gain. Also, p ,' s (for i = l,2) are 
subjected to: 

” P max. ~ Pi ^Pmax ( for 1 = 1,2), (A6) 

for a certain upper bound denoted p max . The respective 
unsteady heat effect Q j associated with these quantities p. 
(for i = 1 ,2 ) are then approximated by: 


Qj =k 


Pi, 


f Po 


+1 -1 


(for i = 1,2) 


(A7) 


In Equation (A7), P Q is a reference pressure comparable to the 
mean pressure drop across the injector. The 
parameters { h . . j , to, K* , p max , and N are chosen to match 

the open- loop forced response characteristics of the 
combustion chamber pressure. 


Since k (in Equations (A4) to (A7)) is small relative to N , it 
can be set to zero, so that Q would be treated as a known input 
at time n. But, the large time-index shift from S N + n to the 
feedback y n is approximated by a 5th-order Pade 
approximation. That is, y n is viewed as the result of a 5 th - 
order Pade-transform (for a time-delay ~ NxT, T = 0.0002 
seconds) applied on the 2x1 vector IVL S v + W . Consequently, 
instead of Equation (Al-b) the following state-space 

relationship between of S N + n and y n should be used with 
(Al-a) to derive a Kalman “estimator” (predictor) for S N + n \ 


Z„.,= UZ„ + FM$ v+n 

y„=EZ„+DMS v+ „+V 


where V is a band-limited random vector; and U, F, E, D are 
the combined state matrices of the pair of Pade- 
transformations (with the respective dimensions 10x10, 10x2, 
2x10, and 2x2). 
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